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A planetesimal moving in the Solar Nebula experiences an aero- 
dynamic drag which causes its orbit to circularize and shrink. 
However, resonant perturbations from a protoplanet interior to 
the planetesimaPs orbit can counteract both the orbital decay and 
the damping of the eccentricity: the planetesimal can be captured 
into an orbital resonance and its eccentricity pumped up to a 
modestly high equilibrium value. Thus, orbital resonances consti- 
tute (partial) barriers to the delivery of planetesimals into the 
feeding zone of the protoplanet. We have established the character- 
istics of the phenomenon of resonance capture by gas drag in the 
circular restricted three-body approximation. We have determined 
the strengths of the equilibrium resonant orbits with respect to 
impulsive velocity perturbations. We conclude that planetesimals 
captured in orbital resonances are quite vulnerable to being dislo- 
cated from these orbits by mutual planetesimal interactions, but 
that the resonances are effective in slowing down the rate of orbital 
decay of planetesimals. Only very small bodies, S100 m, are able 
to reach a ~ lAf 0 protoplanet without being slowed down by 

resonances. CO 1993 Academic Press, Inc 


I. INTRODUCTION 

The currently popular paradigm for the formation of 
the planets (see Wetherill ( 1 989) for a recent review) starts 
with 1- to 10-km-diameter bodies — “planetesimals” — as 
the building blocks. These are presumably the result of 
an early agglomeration process from much smaller solid 
grains present in the Solar Nebula. Large numbers of 
planetesimals are thought to accumulate by means of a 
“runaway” process into a few planetary embryos of 
0(0.01-0. l)A/ 0 in the terrestrial planet region and 
O(1-5)M 0 at the present orbit of Jupiter. At the end of 
the runaway, these planetary embryos have depleted the 
supply of solid material in their “feeding” zones and are 
to be found in relatively isolated orbits. Further growth 

Presented at the Planet Formation program held at the Institute for 
Theoretical Physics of the University of California at Santa Barbara, 
in Fall 1992. 


requires a means of delivering material to the feeding zone 
and/or direct merger of the runaway bodies. 

Most of the recent calculations of the formation of plan- 
etary embryos have neglected the presence of nebular 
gas. In the early stages of the accumulation of the planets, 
solid bodies undoubtedly coexisted with the dominant 
gaseous component of the Solar Nebula. In particular, the 
accumulation of the cores of the jovian planets occurred 
before the dispersal of the nebular gas. Planetesimals 
moving through the gas would experience an aerodynamic 
drag that would damp their random motions (i.e., eccen- 
tricities and inclinations) and cause their orbits to spiral 
into smaller heliocentric distances (Adachi ct al. 1976). 
Naively, then, one would expect that this inward flow of 
solid material would enhance the mass supply rate to the 
accretion zone of a planetary embryo embedded in the 
Solar Nebula. However, this does not necessarily occur 
since planetesimals would also be subject to mutual gravi- 
tational interactions. In particular, resonant gravitational 
perturbations from the protoplanet can counteract the 
effects of gas drag forces, leading to orbital resonance 
locks. A planetesimaPs orbit would then cease to decay, 
and its major axis and eccentricity would attain equilib- 
rium values. This was first pointed out in Weidenschilling 
and Davis (1985). The phenomenon is similar to the 
“mode-locking” phenomena that occur in many other 
natural systems of coupled oscillators. 

The resonance trapping phenomenon acts as a barrier 
to the planetesimals’ inward flow induced by the gas drag. 
However, it is only a partial barrier. In this paper, we 
describe the results of a systematic study of the permeabil- 
ity of this dynamical barrier. In order to elucidate the 
basic phenomenon, we consider the simplest case of the 
circular restricted three-body problem and use a very 
simple prescription for the gas drag force. We have carried 
out extensive numerical simulations to determine the or- 
bital evolution of a planetesimal upon passage through a 
first-order orbital resonance with a protoplanet interior 
to its orbit. We describe the conditions for capture into 
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resonance and the characteristics of the equilibrium reso- 
nant orbit. The numerical results have been compared 
with analytical theory for adiabatic resonance passage 
(Henrard and Lemaitre 1983, Borderies and Goldreich 
1984, Dermott et at. 1988, Malhotra 1988). We find that 
the theory gives reliable predictions for the evolution. In 
a more realistic scenario, planetesimals will suffer gravita- 
tional encounters and inelastic collisions with other plane- 
tesimals. The role and effectiveness of resonance capture 
in the presence of these processes needs to be examined. 
As a crude approximation, we simulate the effect of a 
gravitational encounter or collision by a velocity kick, <5v, 
imparted to the planetesimal. In this manner, we have 
determined the strength and weakness of a resonance 
with respect to planetesimal-planetesimal collisions. We 
demonstrate that in general, a planetesimal may be 
knocked out of an orbital resonance with a relatively small 
dv imparted in the tangential direction, whereas the reso- 
nance is stable against much larger Su’s imparted in the 
radial direction. We estimate that velocity kicks due to 
mutual planetesimal interactions are likely to exceed the 
strengths of resonances; as a result, the resonance traps 
are effective only in slowing down orbital decay, not com- 
pletely halting it. 

The paper is organized as follows. In Section II we 
describe the physical model we assumed. Section III sum- 
marizes relevant analytic estimates obtained from the adi- 
abatic theory of resonance passage in the circular re- 
stricted three body problem. In Section IV we describe 
the results of numerical simulations and compare them 
with the analytic estimates. In Section V we discuss the 
strengths and weaknesses of orbital resonances with re- 
spect to planetesimal encounters. In Section VI we sum- 
marize our conclusions. 


II. THE PHYSICAL MODEL 

We assume that a planetesimal moves in the midplane 
of the Solar Nebula in a heliocentric orbit in the presence 
of nebular gas of constant density and is perturbed by a 
protoplanet which also moves in the midplane in an orbit 
interior to that of the planetesimal. We assume the re- 
stricted three-body formalism and neglect any evolution 
of the protoplanet’s orbit. In heliocentric coordinates, the 
equations of motion for the planetesimal’s position and 
velocity vectors are given by 


r = v, 


v = 


-(1 - ia)^~ n 
r 



f d - 


(i) 


We choose units such that the gravitational constant 
G = 1, and M ( . } + M p = 1 (where M 0 and M p are the 


masses of the Sun and the protoplanet, respectively). The 
parameter p, is the mass of the protoplanet in units of the 
total mass M . ; + A/ p , 


M Mr + M ' 

P 


( 2 ) 


r p is the heliocentric position vector of the protoplanet, 
and f d is the gas drag acceleration; in general it opposes 
the velocity of the planetesimal relative to the gas, but 
depends essentially upon the character of the relative 
motion (Landau and Lifshitz 1987, Adachi et al. 1976). 

We will adopt, for specificity, the following numerical 
values of the physical parameters of the primordial Solar 
Nebula in the vicinity of Jupiter’s orbit. (Although our 
considerations have wider applicability, we focus on these 
conditions at Jupiter’s orbit because planetary embryo 
accumulation here undoubtedly occurred in the presence 
of nebular gas.) These values are used in the subsequent 
sections without further reference. The distance of the 
protoplanet from the Sun is assumed to be a p = 5 AU = 
7.5 x 10 n cm, and its circular orbital velocity is u kep = 
1.3 x 10 6 cm sec" 1 . An ambient nebular gas density 
Pgas - 10“ 10 g cm" 3 and temperature , T = 300 K, are 
assumed. At this temperature, the mean thermal speed 
of gas molecules is c = (3kT/7rm H y f2 « 2 x 10 5 cmsec“ I . 
The molecular viscosity of the gas is then given by v = 
m H /73cr ~ 1 x 10 4 g cm" 1 sec" 1 (where we have used 
<x = 2 x 10“ 1? cm 2 for the collision cross section of 
molecular hydrogen). The gas is assumed to be in near- 
Keplerian rotation. The pressure support in the gas causes 
the gas velocity to be slightly smaller than the local circu- 
lar Keplerian velocity: 1 - (u gas /u kep ) * (l/2)(c7u kep ) 2 * 
5 x 10 The bulk density of the planetesimals will be 
taken to be p = 2 g cm” 3 . 

Taking the relative velocity, //, of a planetesimal with 
respect to the gas to be about 1% of the Keplerian velocity, 
we have u ~ 1.3 x 10 4 cm/sec. (Note that this relative 
motion is subsonic, but just barely so!) For the relative 
motion of a planetesimal of radius s through the gas, the 
Reynolds number, Re = Isp^idv, has a simple mnemonic 
form. 


Therefore, for planetesimals having radii greatly ex- 
ceeding 10 m, Re > 10. At high values of Re, it is an 
experimentally observed fact that the drag acceleration 
is independent of the viscosity and is well represented by 
(Landau and Lifshitz 1987) 

f d = -Kim: u = v - v gas , (4) 

where the coefficient K is given by 
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„ C>s 2 3Cjp gas 

K ~ ~2m~ Pgdi ~ 


(5) 


where m is the mass of the planetesimal, and C d is a 
nondimensional coefficient whose magnitude is known 
from laboratory experiments to be nearly constant at a 
value of about 0.5 for Re in the range 5O0-10\ Thus, for 
example, the numerical value of K is 0.014/AU for a 
planetesimal 100 m in radius. We note that for gas densi- 
ties smaller than assumed above, or smaller sizes of plane- 
tesimal, or smaller relative velocities, the Reynolds num- 
ber could well drop below the range of applicability of 
Eq. (4). If Re ^10, then Stokes 1 law would provide a 
better description. In this paper we confine our attention 
to a drag force in the form of Eq. (4). 

In Eq.(4), u is the velocity of the planetesimal relative 
to the gas velocity, v gas , which we assume is given by 

v gas (r) = (1 - i 7 )v kep (r), (6) 

where 17 ~ 0.005 measures the slight offset (due to the 
partial pressure support) of the gas velocity with respect 
to the local circular keplerian velocity. 

The equations of motion (Eqs. ( 1 )) were integrated using 
a simple modification of a mixed variable symplectic 
scheme (Malhotra 1993), using at least 20 steps per orbit 
of the planetesimal. This numerical method is about an 
order of magnitude faster than conventional integrators 
such as the Bulirsh-Stoer integrator (Press et al. 1986). 
The reliability of the results was confirmed in a large 
number of cases by repeating the integrations using the 
Bulirsh-Stoer integrator. 


III. ANALYTICAL ESTIMATES 

The evolution of a keplerian orbit due to gas drag has 
been analyzed in Adachi et al. (1976). To lowest order in 
the small quantities e and 17, the variation of the orbital 
semimajor axis and eccentricity under the drag law of Eq. 
(4) are given by 


(^)drag (0.97<? + 17)7 


(^)drag “ (0.77e + 7/)e, 

T o 


(7) 


where 


2 x 10 6 (s!\ km) years, and the orbit circularization time 
scale is r e ~ 2 x 10 4 (s/1 km) years. 

In the context of the circular restricted three-body prob- 
lem, in the vicinity of a j + 1: j exterior resonance with 
the protoplanet, the critical resonance angle is defined by 

<t> = U + 1)X - jk p - tu, (9) 

where k and k p are the mean longitudes of the planetesimal 
and protoplanet, respectively, and vj is the planetesimaPs 
longitude of perihelion. The resonant gravitational torque 
from the protoplanet induces a variation of a and e given 
by 


(^)res = -20' + 1 )finaef(a, y)sin </>, 
(Ores = -pnf(aj) sin (t>. 


( 10 ) 


where n = a " 3/2 is the mean motion of the planetesimal; 
a = I la, and/(a, j) is a coefficient of the resonant term in 
the expansion of the protoplanet’s gravitational potential, 
given in terms of Laplace coefficients, by (Brouwer 

and Clemence 1961) 


/(«, j) = | (lj + 1 + a £) //''(a). (11) 

At exact resonance, a = a res = (j/(j + l)) 2/ \ we have 
the asymptotic expression, 

/(a res ,7) = 0.80 j, for j > 1. (12) 

For small rates of orbit evolution due to drag (Eqs. 
(7)), we can draw upon the adiabatic theory of resonance 
passage to determine the conditions for resonance cap- 
ture. A review of this analysis is given in Peale (1986). 
In this theory, the three-body Hamiltonian is averaged to 
obtain a pendulum-like resonance Hamiltonian with one 
degree of freedom (corresponding to the resonance angle, 
</>). The drag-induced evolution is modeled as a drift of the 
particle across the slowly varying phase space trajectories 
determined by the resonance Hamiltonian. As the particle 
energy approaches that of the exact resonant orbit, a 
separatrix appears in the phase space which defines the 
resonance and nonresonant regions. Depending upon the 
initial orbit of the particle and the rate of drag-induced 
drift, the particle may be captured into the resonance 
region or may continue to evolve past the resonant orbit. 1 


T 0 = „ 1 7 7 . (8) 

KUkepM 

For relative velocity u on the order of 1% of the local 
keplerian velocity, the orbital decay time scale is r. d ~ 


1 We note that this theory neglects the damping of the orbital eccentric- 
ity by the drag force. This is a shortcoming that is not of much import 
for the approach to the resonance, but it makes the theory inadequate 
for the evolution subsequent to resonance capture when the damping 
rate becomes comparable to the rate of resonant excitation of the eccen- 
tricity. 
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Qualitatively, the conditions for resonance capture re- 
quire that (i) the planetesimal eccentricity upon approach 
to the resonance be not too high, and (ii) the rate of 
decay of the orbit be sufficiently slow. We quantify these 
conditions below. 

For the first condition, the critical eccentricity is given 
by (cf. Henrard and Lemaitre 1983, Borderies and Gol- 
dreich 1984, Dermott et al. 1988): 


e rr \t = 


~ 2V6 
O' + D 2 


" 1/3 

- 1.58 



La + D 2 J 


1/3 


(13) 


If the second (“adiabatic”) condition is satisfied (and we 
quantify that below), then capture into the resonance is 
certain for e < e cnr For higher eccentricities, the capture 
probability drops off rapidly: 


where e is the initial eccentricity long before the approach 
to the resonance. If the adiabatic condition is only weakly 
satisfied, then capture into resonance can still occur, but 
with probability smaller than that corresponding to the 
adiabatic case (Malhotra and Dermott 1990). Since plane- 
tesimal eccentricities are damped efficiently by nebular 
gas drag, the condition e < e cra is likely to be satisfied 
in many cases. Then, only those planetesimals that can 
overcome the adiabatic threshold may pass through the 
resonance barrier. 

The second condition for resonance capture may be 
quantified by requiring that the change in the semimajor 
axis due to gas drag in a characteristic resonance libration 
period, T t , be smaller than the amplitude (Au)< of resonant 
perturbations of the semimajor axis (Dermott et al. 1988): 



_ 4.0 aKr)(0.91 + V /e) 

^ c; d JO'+i) 


(18) 


such that for a given planetesimal (hence K ), a protoplanet 
of mass (jl > /x th is required to halt the orbital decay of 
the planetesimal by capture in the j + 1 : j resonance. 
C ad is a numerical constant that measures how well the 
inequality (15) needs to be satisfied for resonance capture 
to occur (i.e., C ad is the ratio of the left-hand side to the 
right-hand side). We call it the adiabatic parameter. 

We note that /z (h defined above depends upon the plane- 
tesimaPs orbital eccentricity as it approaches the reso- 
nance. As mentioned previously, if the initial eccentricity 
is too large, resonance capture is unlikely to occur, and 
no other conditions are relevant. However, orbit circular- 
ization by gas drag is quite rapid; therefore, we must 
consider the case of low-eccentricity orbits. If the initial 
eccentricity far away from the resonance is very small, 
then upon approach to the resonance, there is an increas- 
ing eccentricity forced upon the planetesimal. An appro- 
priate value of e to use in the expression for fx tb above is 
the magnitude, e*, of the forced eccentricity when the 
resonance separatrix first appears in the phase space. This 
value is (cf. Malhotra 1988) 

= (2/3) 1/2 ^ cnt , (19) 

where e crit was given in Eq. (13) above. For the parameter 
ranges of interest, we find that the magnitude of e* is 
likely to be comparable to 77 . Therefore, we will assume 
(0.97 + rj/e*) « constant, and absorb its value in the 
numerical factor, C ad = C ' d /(1 + 77 /e); thus, we arrive 
at the following expression for /x th : 


4.0 aK v 
M,h " C ad j(j + !)• 


( 20 ) 


(a)dragT (AaV 
T t and (Aa) f are given by 



9n{nT ( )~ l 
30'+ 1) 


2 

3 


H '2 


M/(« resO> 


We show below that C ad is insensitive to the system pa- 
rameters, and that its numerical value is ^ 3.3. 

In the evolution subsequent to capture into resonance, 
the planetesimal settles into an equilibrium orbit. This is 
.. easily seen by considering the total time rates of change 
ofaand<?: 


From Eqs. (7), (15), and (16), the adiabatic condition is 
found to be: 


aKr) ( 0.97 + <0.25j(j + 1) /x. (17) 


We can define a threshold mass. 


a = -2 (j + \)fjinaef(a,j)s\n </> (0.97^ + 17 )iqa, 


( 21 ) 


e = - finf(a,j)sin <j> (0.77^ 4- r))e. 

T o 


The equilibrium resonant orbit is given by setting a and 
e to zero. This condition yields 
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sin 0 eq 

where corrections of have been neglected. The 

equilibrium eccentricity is independent of the drag param- 
eter, K, and hence of the planetesimal size. For vanish- 
ingly small drag rates, K —> 0, we recover the classical 
result that the fixed points for a first-order resonance are 
at <f) = 0 and <£> = 7r; a linear stability analysis shows 
that the latter is the center of libration in the resonance. 

We can estimate the time scale to establish the equilib- 
rium state as follows. Immediately upon capture into reso- 
nance, we have (a) = 0, where ( ) denote a time average 
over a characteristic resonance libration period. There- 
fore, from Eqs. (21), 


- 1.12 


- 1.21 


T^T. 

J + 1 
7}aK 

j(j + Om’ 


( 22 ) 


(sin <j>) = 


1 (0.97? + r})ig 

TqU + I )pnef(ajY 


(23) 


Inserting this into the equation for ? , we find that the 
secular variation of the eccentricity upon resonance cap- 
ture is given by 



To _ 


(0,97? + 77)77 
O’ + \)e 2 


(0.77? + tj) . 


(24) 


The two terms on the right-hand side represent the reso- 
nant excitation and the gas-drag-induced damping, re- 
spectively. An explicit expression for (?) as a function of 
time is not easy to obtain from this deceptively simple 
equation. However, the time scale for the resonant excita- 
tion of the eccentricity is given by the inverse of the first 
term. By setting ? = ? eq in that expression, we estimate 
the time to establish the equilibrium state to be 


parameter to the system parameters (planet mass, orbital 
elements of the particle, and the parameters of the gas) 
holds only for the first-order resonances that we consider 
here and is a consequence of the form of the drag law. 

The estimates obtained above for the equilibrium reso- 
nant orbit (Eqs. (22)) are equivalent to those found pre- 
viously by Weidenschilling and Davis (1985). (Beware, 
however, of typos in their paper). Those authors argued 
that the condition |sin </> | < I defines a maximum value, 
/C max , for the stability of the equilibrium resonant orbit 
and determines the condition for resonance capture. In 
the notation used here, we get 

«*max - 0 + ^ /Kes J). (27) 

V 

This expression is similar to the adiabatic threshold deter- 
mined by our Eq. (20) above. 2 

In the next section we compare the analytic estimates 
made above with the results of numerical calculations. 
We note here that this analysis is valid only for resonances 
sufficiently well separated that a single resonant term in 
the protoplanet’s gravitational potential dominates the 
perturbations felt by the planetesimal. Sufficiently close 
to the protoplanet, first-order orbital resonances will not 
be well separated, and the perturbations will be a superpo- 
sition of many resonances simultaneously. The response 
of a particle to overlapping resonances is chaotic (Chiri- 
kov 1979). The value, J c , such that all j - f- 1 : j resonances 
with j > j c in the vicinity of the protoplanet are overlapping 
is given by (Wisdom 1980) 

7c - 0.5fx 2n . (28) 

Therefore, we may with some justification consider the 
feeding zone of the protoplanet to extend to the Jc + 1 ; 
7 C resonances interior and exterior to its orbit, i.e., a half- 
width 


Teq ~ (T 0 /e tq ) 


~ 3 x 10 2 (7 + 1) 1/2 



years. 


(25) 


— - 1 .SfJL 217 . 


(29) 


For ju, < 4 x 10 5 , the feeding zone is larger than the 
—3.5 Hill radii that is often assumed. 


r eq is very short compared to the orbital decay time scales. 

A useful exercise at this point is to combine Eq. (20) 
with the expression for sin 0 eq (Eq. (22)) to obtain the 
equilibrium value of </> at the resonance capture threshold: 

sin (/> eq th = -0.30C ud . (26) 

Now, because |sin <j>\ < 1, it follows that C ad s 3.3. It is 
worth pointing out that the insensitivity of the adiabatic 


2 We should point out here a technical difference between our calcula- 
tions and those of Weidenschilling and Davis (1985) in arriving at the 
above condition for resonance capture. These authors implicitly as- 
sumed e — C(tj ) 1 ' 2 > 17 ; this is a reasonable approximation for the equilib- 
rium resonant orbit. However, whether capture occurs is determined by 
conditions at the approach to the resonance. We have argued above 
that, at the approach to the resonance, e is comparable to 77 . Since the 
factor r }le appears only as (0.97 + rjie) in the adiabatic condition, our 
results are not different, but they go to the heart of the mechanism of 
resonance capture. 
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time 

FIG. L An example of capture in 5:4 exterior resonance, (a) The 
equilibrium resonant orbit in the frame rotating with the Sun and the 
perturbing planet: note the characteristic fourfold symmetry (./-fold for 
a j + 1 : / resonance), (b)-(d) The time evolution of the semimajor axis, 
eccentricity, and resonance angle: the unit of time is the orbital period 
of the perturbing planet. The planet mass and drag parameter in this 
example were fx 3 x 10 6 and K = 0.01, respectively. 


IV. NUMERICAL RESULTS 

In Fig. 1 we show the typical orbital evolution of a 
planetesimal upon capture in a first-order resonance with 
a protoplanet which moves on a circular orbit. Initially 
the semimajor axis decays and the eccentricity is damped. 
After the particle is captured in the resonance, the semi- 
major axis stabilizes at the resonant value, a — ((j + 1)/ 
j) m ; the eccentricity and resonance angle approach their 
equilibrium values quite rapidly. We note that long before 
the equilibrium state is attained, the resonance angle may 
be found to be in libration about tt. This libration is dynam- 
ically inconsequential — it arises merely due to the high 
rate of orbital precession induced by a small eccentricity 
in the neighborhood of a first-order resonance, and there 
is no net transfer of angular momentum from the protopla- 
net to the planetesimal during this stage. This is confirmed 
by the observation that the planetesimal eccentricity does 
not get excited and the semimajor axis continues its decay 
unaltered during this stage. 


Figure 2 shows the evolution of a planetesimaFs eccen- 
tricity upon passage through a 5 : 4 orbital resonance for 
several different protoplanet masses. It may be seen that 
for cases where /x is greater than a threshold value, /x th , 
the eccentricity evolution is independent of the value of 
/ x ; the planetesimal is captured in the resonance and set- 
tles into an equilibrium state. For values of /x smaller than 
resonance capture does not occur; the eccentricity is 
briefly excited to a modest value but is eventually damped 
away. The magnitude of this excitation decreases rapidly 
with the mass of the protoplanet. 

In Figures 3, 4, and 5 we have summarized the results 
of an extensive series of numerical simulations of the 
resonance capture phenomenon. We simulated the evolu- 
tion for two different values of the drag parameter: 
tf p K = 0.001 and a p K = 0.01, and a wide range of pro- 
toplanet mass, /x. For reference, these two values of K 
correspond to approximately 7.5- and 0.75-km-radius 
planetesimals, respectively. 

• Figure 3 shows the minimum protoplanet mass neces- 
sary for a planetesimal with a small initial eccentricity 
( e < c crjt ) to be captured into resonance; the dashed line 
indicates the adiabatic threshold defined by Eq. (20). 

• Figure 4 shows the equilibrium eccentricity as a func- 
tion of the integer, j\ the numerical experiments confirm 
its insensitivity to the protoplanet mass and the planetesi- 
mal size. 

• Figure 5 shows the numerically determined equilib- 
rium value of the resonance angle at the capture threshold. 
There appears to be a weak dependence on the particle 
size and the integer ./, such that the center of libration, 



FIG. 2. The time evolution of a planetesimaFs eccentricity upon 
passage through a 5 : 4 exterior resonance. Four cases, each correspond- 
ing to a different perturbing planet mass, /x, are shown. For /i less than 
a certain threshold value, the planetesimal passes through the resonance 
without being captured. 
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J 

FIG. 3. The threshold planet mass for capture into j + 1: j exterior 
resonance. The points are from numerical experiments, the dashed line 
represents the analytical estimate of Eq. (20) with numerical coefficient 
C ad = 3.0 for the case K = 0.01 , and C ad = 2.5 for the case K = 0.001 . 



j 

FIG. 4. The equilibrium eccentricity in j + 1: j exterior resonances, 
as a function of the integer,,;'. The points are from numerical experiments, 
the dashed line represents the analytical estimate of Eq. (22). e cq is 
insensitive to the planetesimal size and to the mass of the protoplanet. 


<£ eq , for smaller planetesimais at distant resonances is 
removed farther from the “classical” value of 7 r. 

These results show that in the restricted three-body for- 
malism, the analytical estimates describe the dynamical 
evolution quite well. 

However, for parameters close to the resonance cap- 
ture threshold, the dynamics is very complicated and the 
analysis of Section III can be relied upon only broadly, 
but not in the utmost details. In particular, we find that 
the analytical estimate of <? eq is quite robust but sin </> eq is 
not quite so. For example, Fig. 5 shows that at the capture 
threshold, the equilibrium values of the resonance angle 
hover close to 37 t/ 2; however, the resolution of our numer- 
ical experiments was too poor to make a meaningful inter- 
pretation of the deviations from 37 t 72. Recent work by 
Kary et aL( 1993) discusses this issue in greater detail. 


V. STRENGTHS OF RESONANCES 

In this section, we consider the strength of a j 4- 1 : 
j exterior resonance against close encounters between 
planetesimais. As a crude approximation, we can model 
the effect of a planetesimal-planetesimal encounter by an 
instantaneous kick in velocity, 8v, imparted to each body 
which, in effect, reinitializes the orbital elements. As long 
as the new orbital elements remain within the resonance 
width, the planetesimal will remain trapped in the reso- 


nance. The resonance width varies sensitively with orien- 
tation in phase space. Therefore, we expect the strength 
of the resonance to be sensitive to the direction of the 
velocity change, 8v. 

2tt 
6 


5 

cr 
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7 T 
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j 

FIG. 5. The numerically determined equilibrium value of the reso- 
nance angle at the resonance capture threshold. 
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FIG. 6a. Time evolution of the semimajor axis and eccentricity of 
a set of 20 planetesimals {K - 0.01 , n — 3 x 10 \ (Unit of time is the 
orbital period of the planet.) All the planetesimals are at first captured 
in the 5:4 resonance. At time = 6103.7, each planetesimal suffers a 
velocity kick of magnitude 8v = 0.005 but different orientation. In the 
subsequent evolution, some remain the 5 : 4 resonance, others are kicked 
out of the 5 :4 but get captured in other resonances closer to the pro- 
toplanet. 



t 

Pianetesima! 


FIG. 6b. The strength of the 5 : 4 resonance with a 10 M® protoplanet. 
Shown is a summary of the results from numerical experiments of the 
type indicated in Fig. 6a: the shaded regions indicate the magnitude and 
orientation of an impulsive velocity kick that does not result in the 
planetesimal being ejected out of the resonance. The position of the 
particle at the time of the velocity kick is indicated by •. It is apparent 
that the resonance is far stronger against radial 8v’s than against tangen- 
tial Sv's. 


We can estimate as follows the maximum dv that a 
planetesimal captured in a resonance can suffer without 
being knocked out of the resonance. A dv which is along 
the velocity vector of the planetesimal will produce a 
change in the semimajor axis, da = 2 (dv/v)a. The maxi- 
mum half-width of the resonance along a is simply (Aa)* 
(Eq. (16)) evaluated for the equilibrium resonant orbit 
(i.e., at e = f eq , Eq. (22)). Thus, we obtain the following 
estimate: 


dv 

v 


- 1.1 


'res 



tl V2 V ]IA 


0 .2?/ ,/ V /2 . 


(30) 


We have verified this estimate by means of numerical 
experiments. An example is shown in Figs. 6a and 6b. 
Figure 6a shows the evolution of a set of 20 planetesimals 
(with the same initial conditions) into a 5 : 4 exterior reso- 
nance with a 10A/ e protoplanet. At a specified time into 
the evolution, each particle was imparted a velocity kick 
of magnitude dv = 5 x 10“ 3 (but different orientation). 
It may be seen that some fraction of the particles remain 
in the resonance, or are knocked into higher orbits which 
are then recaptured into the 5 : 4 resonance; the rest are 
knocked out of the resonance and their orbits decay in- 
ward. (In the example shown, all of these eventually get 
captured into other resonances interior to the 5 :4 reso- 
nance.) Figure 6b summarizes the results of many such 
numerical experiments. It shows the strength of the reso- 
nance in terms of the magnitude and orientation of the 
velocity kicks necessary to disrupt the resonance: there 
is clearly a high degree of anisotropy. The theoretical 
estimate of Eq. (30) yields dv/v — 0.002 for the width 
along the tangential direction; this is in very good 
agreement with the results shown in Fig. 6b. 

We have performed similar experiments for several dif- 
ferent resonances (and with differing masses of the pro- 
toplanet) and have found that Eq. (30) provides an excel- 
lent estimate of the resonance strength with respect to 
velocity kicks along a tangent to the planetesimafs orbit. 
Along the radial direction, the resonance strength is 
greater by factors of ~4 or more. 


VI. CONCLUSIONS 

The conditions for resonance capture imply a high effi- 
ciency of the phenomenon for low-eccentricity planetesi- 
mals. In fact, it is interesting to note that a 1 Earth mass 
protoplanet provides resonant torques of sufficient magni- 
tude to halt the orbital decay of planetesimals ~3 km or 
larger in radius even in distant resonances such as the 
3:2. Because resonances closer to the protoplanet pro- 
vide larger resonant torques, a planetesimal that is too 
small to get captured in distant (small j) resonances may 
eventually be captured in a resonance closer in. 
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Given the high likelihood of resonance capture, the 
stability of the resonance traps with respect to mutual 
planetesimal interactions is an important consideration 
for planet formation theories. In the previous section, 
we provided a quantitative measure of this stability. A 
detailed discussion of the implications is beyond the scope 
of the present work (see Torbett et al. 1982, Weidenschil- 
ling and Davis 1985 for qualitative discussions on this 
point), but as an illustration, let us consider the specific 
case of a 1 A/ 0 protoplanet at 5 AU and a swarm of 
10-km-radius planetesimals. From Eq. (30) we have 
(8u) res ~ j i/ 4 6.5 m/sec. Now, in the planetesimal swarm, 
random velocities of encounter are on the order of the 
escape velocity from the surface, which has a similar 
magnitude, —10 m/sec. Furthermore, the relative velocity 
of encounter between a resonant and a nonresonant plane- 
tesimal would be even larger, — e eq * u kep =» (j + l) _l/2 * 
1 km/sec, easily exceeding the strength of the resonance. 
Thus, it is apparent that a planetesimal trapped in a reso- 
nance is quite vulnerable to being knocked out of the 
resonant orbit. However, there are a couple of points to 
keep in mind: (i) unless it is fragmented, a planetesimal 
knocked out of one resonance has a high probability of 
being recaptured in the same or neighboring resonance; 
(ii) if fragmentation does occur, the bigger fragments will 
tend to be imparted smaller 8t/s and so will tend to stay 
in resonance. We conclude that mutual planetesimal inter- 
actions make the resonance barriers permeable, but the 
resonances may be effective in slowing down the orbital 
decay. 

The above considerations apply to those planetesimals 
that do get trapped in resonances. Planetesimals that are 
sufficiently small, 5 < that their rate of orbital decay 
exceeds the adiabatic threshold for all j + 1: j exterior 
resonances with j > j c (Eq. (28)) will be delivered to 
the feeding zone of the protoplanet without being slowed 
down by the resonance traps. By combining Eqs. (5), (20) 
(with C ad = 3.3), and (28), we find the critical size, s Q to 
be 



It is appropriate to recall here that the gas drag law 
we have adopted (Eqs. (4) and (5)) holds with C d ~ con- 
stant ~ 0.5 for high Reynolds numbers. From Eq. (3), we 
see that this assumption may be in question for planetesi- 
mals smaller than the s () determined above. While the 
form of the drag force is expected to hold for Re as small 
as ©(10), the numerical value of C d increases steeply with 
decreasing Re (Landau and Lifshitz 1987). Therefore, we 


have indicated explicitly in Eq. (31) the dependence of ,v 0 
on C d . In any case, it may be concluded that the critical 
size s 0 is on the order of 100 m for a 1 protoplanet. 

In summary, we have discussed the trapping of plane- 
tesimals in orbital resonances with a protoplanetary em- 
bryo in the Solar Nebula. Several idealizations were nec- 
essary for elucidating this fundamental dynamical 
process. Our model included only a single, isolated plane- 
tary embryo, assumed to be in a circular orbit, and treated 
the planetesimals as “test particles"; a constant and uni- 
form nebular gas density was assumed; the response of 
the gas to the gravity of the protoplanet and planetesimal 
was neglected. Within this idealized model, we have es- 
tablished the conditions for resonance capture, the stabil- 
ity of the resonances against mutual planetesimal pertur- 
bations, and the implications for the delivery of solid 
bodies to the planetary embryo. However, the idealiza- 
tions are sufficiently severe that much further work is 
necessary before the influence of the resonance trapping 
phenomenon on the late stages of growth of planetary 
embryos in the Solar Nebula can be fully evaluated. 
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